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The production of electron-positron pairs with the capture of the electron in an atomic orbital is 
investigated for the conditions of the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron 
Collider LHC). Dirac wave functions for the leptons are used, taking corrections to orders of Za 
into account. The dependence on the transverse momentum transfer is studied and the accuracy of 
, the equivalent photon approximation is discussed as a function of the nuclear charge. 
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■ I. INTRODUCTION 

(N : 

, Peripheral collisions with relativistic heavy ions are a source of several interesting physical processes and have 
^ , attracted a great number of theoretical and experimental work (see, e.g., and references therein). One of the 
processes of interest is pair production by the strong electromagnetic field of large Z nuclei in relativistic heavy ion 

' colliders. In the case of production of free pairs the theoretical interest resides in the non-perturbative character of the 

00 ' . . 

' process. Different techniques have been used to study this problem theoretically, but we will not discuss them here. 

Another case of general interest is the production of pairs in which the electron is captured in a bound state of one of the 
colliding nuclei (see also ref. and references therein) . The interest here arises from its consequence for luminosity 
, changing of the relativistic ion beams in colliders like RHIC (Relativistic Heavy Ion Collider) at Brookhaven, and the 
LHC (Large Hadron Collider) at CERN (for a recent analysis, see Q). A more striking application of this process 
^ was the recent production of antihydrogen atoms using relativistic antiproton beams In that case the positron is 
^ produced and captured in an orbit of an antiproton. 
. ^ Early calculations of pair production with capture used perturbation theory , the only difference being the 

way the distortion effects for the positron wavefunction were taken into account. Evidently the calculation is best 

■ performed using the frame of reference of the nucleus where the electron is captured. Many other calculations have 
been performed Some of them used non-perturbative approaches, e.g., coupled-channels calculations. However, 
in contrast to the production of free pairs, pair-production with atomic capture of the electron is well described in 
perturbation theory. The large coupling constant for large nuclei (Za ^ 1) does not matter here, since the scaled 
matrix element divided by Za, i.e. M / {Za), is much smaller in this case than it is for the production of free pairs. In 
other words, the overlap between the positron and the electron wavefunctions is much less than that for the production 
of free pairs. The first terms of the perturbation series are small enough to neglect the inclusion of higher order terms 

. Indeed it was explicitly shown in ref. that non-perturbative calculations modify the cross sections at most by 
1%. 

The theoretical difficulty here is to properly account for the distortion of the positron wavefunction in the field 
of the nucleus where the electron is captured. In this article we study this property more closely, and compare our 
results with the calculations of ref. which were based on the Sommerfeld-Maue wavefunctions for the positron 
. We also investigate the accuracy and limitations of the equivalent photon approximation (EPA) , what will gives 
us a clear understanding of the effects of distortion. We have closely used the technique of ref. il^ for the study of 
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anti- hydrogen production cross sections at LEAR/CERN. We show that a commonly used approximation obtained 
in ref. |^ yields results up to a factor 5 too small for large ion charges. 

In section 2 we deduce the cross sections for pair production with K-shell capture and obtain the cross sections in 
terms of longitudinal and transverse components. In section 3 we compare the calculation with that of pair production 
with capture by real photons. This allows us to disclose the role of distortion effects. We also show that the transverse 
part of the virtual photon cross section is much larger than the longitudinal one. In section 4 we calculate total cross 
sections and compare them to the results obtained in ref. to the equivalent photon approximation and to recent 
experimental data . Our conclusions are given in section 5. 



II. PAIR PRODUCTION WITH K-SHELL CAPTURE 

In the semiclassical approach we assume that the relativistic ions move along straight lines. Although the laboratory 
system is a more symmetric frame of reference, the calculation is best accomplished in the frame of reference of one 
of the ions - the one in which the electron is captured. The reason is that the effects of distortion of the positron 
wavefucntion are very important, and cannot be neglected. This distortion effect is most naturally taken into account 
in the rest frame of the nucleus. 

The relationship between the Lorentz contraction factor associated with the relative velocity between the colliding 
nuclei, 7, and the collider energy per nucleon in GeV , E/A, is given by 7 = 2(1.0735 x E/A)^, for E/A 3> 1. For 
example, at the RHIC, E/A = 100, and 7 = 2.3x lO"*, while at the LHC, E/A ^ 3x 10^ and 7 ^ 2.07x 10^. According 
to refs. the probability amplitude for pair production with capture for a collision with impact parameter b is 

given by 

where qt is a transverse momentum vector with two components, Q = (qt, w/v) is a 3-component momentum vector 
with the third component equal to oj/v, and w is the total energy used in the pair production. It is important to 
keep the relative velocity between the colliding nuclei, v, in the right places for the moment (although v = c = 1) 
, since important combinations of 1 and v will lead to Lorentz factors 7 = (l — v^) in subsequent steps of the 
calculation. In the case of K-shell capture w = £ + m — I^^, is the ionization energy of the K-shell electron, and e 
is the positron total energy e = y/p^ (from now on we will use natural units, so that fi, = c = to = 1). The form 
factor F(Q) is given in terms of the lepton transition current jp(r) — e^'(+'7p5'(~) as 

F{Q) ^leK^j d\ *W(r) e*'^'"^ 7^ *(-)(r) , (2) 

where = (0,K), with K — (qi/oj, l/v7^). For details of derivation of eqs. (|l]^), see refs. [§,0. The cross section 
is obtained by integrating the square of the expression (|^) over all possible impact parameters: 

spins spins 

For the positron wavefunction we use a plane wave and a correction term to account for the distortion due to the 
nucleus charge. The correction term is considered to be proportional to Ze^. The wavefunction is then given by 

= N+ [v e'P " + , (4) 

where iV+ is a normalization constant. 
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Only the Fourier transform of will enter the calculation. This Fourier transform can be deduced directly from 
the Dirac equation for the positron in the presence of a Coulomb field of a nucleus: 



tp.r 



(5) 



Applying on both sides of this equation the operator (j'^s + i'j' —m^ we get 



,ip.r 



47r 



Multiplying by e and integrating over cfir we get 

(p2 - q^) ^'^ = -Ze^ - 7^.(q - p)] (^M , 

where we have used the identity 7" ^7*'e — +rnj (7"v) = 0. Thus, 



(6) 



(7) 



_ 270^+ 7^.(q-p) 



(q- p)2(q2 - p2) 
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(8) 



We will need this result later to calculate F{Q). The normalization constant 7V_|_, which accounts for the distortion 
of the wavefunction to all orders is given by M 



N+ = exp [7ra+/2] F (1 + ia+) , a+ = Ze'^/v+ , 

where v+ is the positron velocity in the frame of reference of the nucleus where the electron is capture. 
For the electron we use the distorted wavefunction described in ref. Ml, i.e.. 



1-1/7^.^ 



u* 



non—T 1 



where 



1 



3/2 



—= cxp(-Zr/a/f) 



with aji — l/e^ being the hydrogen Bohr radius. 

Using the positron and the electron wavefunctions as given by eqs. (0) and (foh, we get 



F(Q) = ieN+ j d\ {v e-'P " + (7^^^) e' 



Q.r 



1 - ^ 7° T"-"^ \ n-fnon-ri^) 



+ 4"e'Q -- (7^i^^)u*„o„_,(r)l , 



U*„on-r(r) 



(9) 



(10) 



(11) 



where in the last equation we neglected terms of highest orders in Za. 
Integrating by parts yields 



1 



1 + ^ I'^-iQ - P) ) U (*no„-,-)Q_p 



1 / z 



3/2 



(12) 



3 /2 

Since the corrections are to first order in Za, we have replaced ^'„o„_r(r) by its constant value {Z/an)' {I/^/tF) 
in the last term of eq. ( p^ . If we would do the same in the first term it would be identically zero (except for Q = p) . 
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We need to keep the corrections in the wavefunctions at least to first order in Za\ for e ^ 1 these corrections yield a 
term of the same order as the plane-wave to the total cross section, as we shall see later. 
Inserting these results in eq. |[ and using [exp (— 'r/a)]q„p = Stt/o (Q — p)*, we get, 

(13) 



aH 



(Q-pr 



where 



with 



(14) 



(Q P)^ Q2 P^' 
Q p 

2(Q2-p2) ■ 



Q-p 

2(Q-P)' 



and 



(15) 



Inserting (|T^ into (|^), with l-ZV+j^ — 27ra+/ [exp (27ra+) — 1], performing the sum over spins with standard trace 
techniques, and using p^dp — epde {p= |p|), we get 

da 64 Z^e^p f d'^qt 1 



dedn 7raf^v+ [exp {2TrZe^/v+) - I] J (Q2 - cj2)2 ^.q _ 
X i (e - 1) [(b - df K^ + A (K.b) (K.d) 



+ (e + 1) a^K'^ - 2a [2 (p.K) (b.K) + p.(b - d)K'^] 



(16) 



In this relation we have used e_ 1, neglecting the electronic binding energy in the numerator of the expression 
inside the integral. 

Using the definition of /^^ , it is illustrative to separate the integrand of the above equation in terms of longitudinal 
and transverse components: 



dedO. TTU! 

where 
dcTj^T 



dn 



(Q,^) 



(Q2-u;2)2 



6„4 



d^ 



dn 



dn 



Z'=e 



(e-l) 



dcT-v* 



dn 



(Q,^) 



\jUJY+ [exp (27rZe2/v+) - 1] (Q - p) 
+ (e + l)a2 + 2a[ ('^--P-^y^ -p.(b-d)1|, 

(Q-p) J J 

M-Kpuj Z^e^ 1 



(b-d) 



{qt - p.et)' 



(Q-p)^(Q2-p2) 



a\qh^^+ [exp (2^Ze2/v+) - 1] (Q - p) 
+ (£ + l)a2 + 2a 



(e-1) 



(b-d)2- 



(w/v-p^) 



(Q-p)' (Q2-p2) 



2 P(o d) 



(Q,c^) 



(Q-P)^ 

1287rp 



(w/v-p^) 



aH«t7'v+ [exp (27rZe2/v+) - 1] (Q - p)^ 



X < a 



(p.et) +p; 



(gt - p.et) 

(tj/v-p^) 



(e-1) 



(gt - p.ef) 
(Q2-p2); ' 



(17) 



(18) 



(19) 



(20) 
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In these equations e = qt/gt is a unit vector in the transverse direction. The cross sections cr^*T, o'7*l, and a~f*LT 
are interpreted as the cross sections for pair production with capture by virtual transverse and longitudinal photons 
and a interference term, respectively. Note that the transverse and longitudinal directions are with respect to the beam 
axis, not with respect to the photon momentum, as is usually meant by this term. Only for 7 3> 1, this definition 
agrees with the definition of transverse and longitudinal virtual photons. However, this separation is very useful, as 
we will see next. 

Now we turn to the relation to pair-production by real photons and the relation to the equivalent photon approxi- 
mation. 



III. PAIR PRODUCTION WITH CAPTURE BY REAL PHOTONS: ROLE OF DISTORTION EFFECTS 



The cross section for pair production with capture real photons is given by 



2TT\Vf,\'Siu-e-l) 



(27r) 



3 ' 



where 



J dV^'(+)(r) e"^ " (e.7^)^'(-)(r) 



(21) 



(22) 



with it = z Lu, where z is the unit vector along the photon incident direction, and e is photon polarization unit vector. 
Using the positron and the electron wavefunctions from the previous section, and performing similar steps, we get 



n 87r\/2 vSu 
2 TV, 



where 



B = a^i^.H) + (t-^)7° ( T'-b.) + ( 7'.d.)7°( 7^-^) 



(23) 



(24) 



where a„, b^, and d„ are the quantities defined in eqs. |l5| but with Q replaced hy it = z oj. Inserting these results 
in eq. |l] and summing over spins we get 



Z^e^ 



1 



da^ . , GAire^p 
In ~ a%LJv+ [exp (27rZe2/v+) - 1] {it-pf 



ie-l) 



(b.-d.)^ 

("^— p) (k^ — p^) 



+ {e + l)al + 2a^ 



iP-ef 



p.(bK d^) 



(25) 



We notice that the above equation can also be obtained from eq. ^ in the limit qt ^ and v — > c. In this limit 
a~f-L — * Oj and (J^y-LT — * 0, and a-^-T becomes the cross section for the production of pairs, with capture, by real 
photons. 

Integrating eq. |2^ over the azimuthal angle we get (without any further approximations!) 

^6g2 g2p3 2) sin^ e 



da. 



de ^'^^ [exp (2^Ze2/v+) - 1] a]jLu^^+ {e-pcosOf 



e — p cos ( 



uj{e + 2) 



(26) 



In the frame of reference of the nucleus, the angular distribution of the positrons is forward peaked, along the 
incidence of the photon. The higher the positron energy is, the more forward peaked the distribution becomes. In a 
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collider, a Lorentz transformation of these results to the laboratory frame implies that all positrons are seen along 
the beam direction, within an opening angle of order of I/7 <C 1. 
Integrating eq. over 9 we get 



P 



1% v+ [exp (27rZe2/v+) - 1] tj* 



P 



■ In (e + p) 



(27) 



The relevance of the Coulomb distortion corrections can be better understood by using the high energy limit, e 3> 1. 
From eq. ( p7|) we get 



1 



v+ [exp (27rZe2/v+) — 1] £ 



(28) 



If we had used plane-waves for the electron and hydrogenic wave function for the positron in our calculation, without 
the corrections to order Za, the cross section for the pair production with capture by real photons would be given by 



(c^) 



= 47r- 



PWA 



(e — pcos( 



the integral of which being 



In the high energy limit: 



SvrZ^e^ p 



''H 



4 (3£^+/) 



5^2 



PWA 



32TrZ-'e 
3< 



(29) 



(30) 



(31) 



Using the approximation Za ^ 1 in eq. 



we get 



(32) 



''H 



which is different than eq. (31) by a factor 8/3. This is important to show how the corrections of order Za influence 
the result of the calculation. Not only they modify the results by inclusion of the correct normalization of the 
distorted positron wavefunction, iV+, but they also yield terms to the matrix elements for photo-production of the 
same magnitude as the terms of lowest order. The origin of the difference between the two results are the small 
distances which enter in the calculation of the matrix elements of eq. (p2|). These corrections are essential to account 
for their effects properly, and are enough to account for a good description of pair production with capture, as we will 
see in the next section. It is worhtwhile to mention that this equation confirms the findings of Sauter ||2^ for the 
annihilation of a positron with an electron in the K-shell of an atom. Using the detailed balance theorem that process 
can be related to the gamma production of electron-positron pairs with electron capture and the above equation is 
reproduced. 

All positron variables in the above equations are expressed in the frame of reference of the ion where the electron 
is captured. Now the (Q^ — o;^)^ denominator of the first term inside the integral |l^ shows that qt = uj/v^ is the 
relevant order of magnitude of qt in the integration. Thus, in all three virtual photon cross sections, 18, 1^, and ^ we 
have IQI = lu/v = lu ^ lu/yj. The additional terms seem also to be independent from 7, what would imply that no 
extra 7 factor arise in these expressions. Thus, for large 7, a-y'T — 7^cr7*L = Jcr^y'LT- This is confirmed by numerical 
calculations. Thus, for RHIC and LHC energies we can safely use 



d(T 1 

Z e 

TTLU 



2„2 



dedVl 



Qt 



(Q2 



da. 



7-T 



dn 



(33) 
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The above discussion might seem convincing. But, it is still reasonable to check this result within a solvable 
model. This can be achieved by using plane-waves for the positron and hydrogenic waves for the electron, without 
the correction terms to order Za. In this case, we get 



(Q-p) 



which yields the cross section 



da 



^ ry2 2 

— Z e 



de(Kl 

where the interference term is exactly zero, and 



(Q2-C^2)2 



da. 



7-T 



d^ 



da. 



7-L 



dVL 



(Q,w) 



dn 



(Q,^) 



(Q - P)' 



and -^^(Q'^) 



2 

7^v 



(Q-pr 



The integral over in the expressions above can be done analytically yielding 



a~f*L {Q,u!) 



USTiZ^e^p (p2 + 3Q2) (3p2 + Q^) / 1 



3a 



H 



2\6 



(Q2 - p2) 



(34) 



(35) 



(36) 



(37) 



Using this result, the integral over in (|35| ) can also be performed analytically but resulting in too long expressions 
to be transcribed here. We also note that the denominator (Q — p)^ in eq. ( [l3| ) implies that |p| ^ |Q| = w/v and 
that the largest contribution to the integrals arise from the region with (Q — p)^ = Thus, in the plane wave 

approximation we really expect that cr-yT is by some power of 7 larger than the (t^»l cross section. 



IV. TOTAL CROSS SECTIONS: COMPARISON WITH THE EQUIVALENT PHOTON 

APPROXIMATION 



The energy spectrum of the positrons can be obtained by integrating ^ over angles. The calculation can be done 
analytically if one neglects terms of order 0(1/7). The expressions obtained were checked against the numerical 
results and we found that for RHIC energies the difference between the two results are only visible (of order of 1%) 
for positron energies e <C 1. For positron energies e 1 the differences are imperceptible. One gets 



da 

de 



327rp 



8^6 



[ CD 

z|jwv+ [exp(27rZe2/v+) ^ \"b6" 



197 
"60" 



3 D 

5 55 



\C_ 



3 1 
20 54 



where 



r 



B 



and D = 3p^ + 



(38) 



(39) 



In ref. Q analytical formulas were also derived for the differential and total cross section of pair production with 
K-shell electron capture. The cross sections were obtained by using Sommerfeld-Maue wave functions for the positron. 
Their result is 



da 
"a 



a%uj^v+ [exp (27rZe2/v+) - 1] 



2e 

y 



P 



In (e + p) 



In (2:^ 



OJ 



(40) 



where d = 0.68108... is related to the Euler's number. 

It is also worthwhile to compare the calculations with the equivalent photon approximation. The function 
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d {ql) - n (7, Q) a^.T {qt = 0, c^) = - {lo) / d (g^) n (7, Q) (42) 



n (7, Q) = 3 Z\^- ?i -2 (41) 

is called by "equivalent photon number", and is strongly peaked at qt — i^/^y, which is very small for 7^1. The 
equivalent photon approximation states that 

da /■^i"'"' , 1 1 /""Jr" 

- n (7, Q) cr^.T (qt = 0, w) = - 
10 ^ ^ JO 

where in the last equality one approximates ct^'t {qt = 0,^^) ~ a-y (uj) . The problem with the approximation is that 
the integral in eq. ^ diverges logarithmically. The approximation is only valid if we include a cutoff parameter qf^'^^, 
determined by the value of qt at which (t^»t (Q,i^) drops to zero. For RHIC and LHC energies, we verified that the 
transverse virtual photon cross section drops to zero at qt — e. Another hint for obtaining the appropriate cutoff value 
of qt is given in figure 2 where we plot the energy spectrum of the positron, da /de, for RHIC energies. We observe 
that the energy spectrum peaks at e ~ 2 — 3. We thus conclude that an appropriate value of the cutoff parameter is 
^max ^ 2. Inserting this value in eq. ^ we get the same result as|4^, except that in the last logarithm of that formula 
5 is replaced by 0.6065..., an irrelevant change, in view of the large value of 7. 

In figure 2 we compare our result, eq. ^ (solid line) with the approximation ^ (dashed line). We observe a good 
agreement, except at the low positron energies. This is an important result. The Sommerfeld-Maue wave functions 
(see fl^ for a complete discussion about these wavefunctions) reproduce the exact Coulomb-Dirac wavefunction in 
the region relevant for the pair-production process. They account for the distortion effects in all orders. The result 
of figure 2 shows that, for pair production with K-shell capture, the approximation ^, with ^ leads to very close 
agreement with that using Sommerfeld-Maue wavefunctions. That is, the most important part of the distortion is 
already taken into account by expanding the wavefunction to first-order in Za and by using the correct normalization 
for the distorted wave ^. We also note from figure 2 that the two results agree at large positron energies, e ^ 1. That 
is expected, since distortion effects are not important at these energies and the outgoing positron is well described by 
plane- waves. 

The total cross section can be obtained from eqs |8| and ^ by an integration over the positron energy. The 
integrations are performed numerically. In ref. 1^ an approximate analytical result was obtained. Their result is 

^^""-Z^e' ^ ^ [In (7<5/2) - 5/3] . (43) 



lQa% exp{27rZe2} - 1 

The term [] ^^ is the main efffect of the distortion of the positron wavefunction. It arises through the normalization 
of the continuum wavefunctions which accounts for the reduction of the magnitude of the positron wavefunction near 
the nucleus where the electron is localized (bound). Thus, the greater the Z, the less these wavefunctions overlap. 

In figure 3 we compare the three results for RHIC energies, as a function of the atomic charge Z. The numerical 
integration of eqs. ^ and ^ yields essentially the same result, differing at most within 2% from each other. They 
are displayed by the solid curve. The dashed curve is the result of It agrees with the integrations of eqs. ^ and 
1^ for low Z (Z < 15) but disa grees by a factor 5 for large Z (Z ~ 80 — 90). The reason for this discrepancy is that 
eq. ^ was deduced in ref. ||] by probably neglecting the strong energy dependence (for large Z) of the exponential 
in the denominator of |3^ which is a function of the positron velocity, v+ . 

We note that the cross section |8| scales approximately with Z'^ for small values of Z, as also implied by eq. 
But, for large Z's there is not such a simple dependence. In our case, it is directly obtained from integration of eq. 
|38| . The charge dependence of the capture cross sections have also been studied in refs. [|l2|,0 with slightly different 
results. For asymmetric systems, e.g., a projectile of charge Zp incident on a target of charge Zt , the equation ^ 
should be changed by replacing Z^ in the numerator by Z^Z^ and Z in the exponential term of the denominator by 
Zx- 
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Equation |8| (and also ^) predicts a dependence of the cross section in the form cr = Alnjc + B, where A and 
B are coefficients depending on the system and bombarding energy This dependence was used in the analysis 
of the experiment in ref. |16|. In recent calculations, attention was given to the correct teatment of the distortion 
effects in the positron wavefunction In figure 4 we show the recent experimental data of ref. compared to 
eq. ^ (dashed line) and to eq. The data was obtained for P6^^+ at 33 TeV impinging on several targets. The 
comparison with the experimental data is not fair since atomic screening was not taken into account here. When 
screening is present the cross sections will always be smaller at least by a factor 2-4 |^ . 

Using eq. |3^ for Au + Au beams at RHIC (100 GeV/nucleon) and Pb + Pb beams at LHC (3 TeV/nucleon), we get 
a — 229 b and a — 562 b, respectively. Other results can be easily obtained from a simple integration of eq. |3^. In 
ref. a numerical calculation was done making use of full Dirac wavef unctions, yielding 114 b for Au + Au beams 
at RHIC. This is about 2 times larger than that obtaining using eq. ^ from ref. ||^ and about 2 times smaller than 
our eq. We assume that the reason for these diferences are that our wavefunctions are correct only to first order 
in Za. It might be that the inclusion of higher-order terms will change our results. However, the nice feature of our 
approach is that it provides an analytical description of the process and a direct connection with the photo-production 
cross sections. 



V. CONCLUSIONS 

We have calculated the cross sections for pair production with K-shell capture in peripheral collisions of relativistic 
heavy ions. The cross section can be separated into longitudinal and transverse components, with respect to the 
beam direction. For high energies, 7 3> 1, this corresponds to the usual meaning of longitudinal and transverse 
components of the photon. A very transparent and simple formulation is obtained using the lowest order corrections 
to the positron and electron wavefunctions. At ultra-relativistic energies, corresponding to the energies of relativistic 
heavy ion colliders (RHIC and LHC), the contribution of the longitudinal virtual photons to the total cross section 
vanishes. The remaining terms of the cross section can be factorized in terms of a virtual photon spectra and the 
cross section induced by real photons. However, the factorization depends on a cutoff parameter, which is not well 
defined. But, for 7 3> 1 a momentum cutoff parameter qt ^ 1 reproduces extremely well the exact value of the cross 
section. 

We note that the electron can be captured in any atomic orbit around one of the ions, the capture to the K- 
shell being the largest contribution. In ref. [|| it was shown that the contribution of capture to all other shells will 
contribute by an increase of 20% of the presently calculated cross section. Depending on the relevance of this process 
to peripheral collisions with relativistic heavy ions and on the accuracy attained by the experiments, more theoretical 
studies on the capture to higher orbits should be performed. At the present stage, it is important to notice that 
eq. ^ from ref. 1^, is not a good approximation for large Z ions. The calculations developed in the present article 
were important to pinpoint the origin of this discrepancy and the validity of the equivalent photon approximation 
for 7^1. It would be nice to show in a future study if the distortion corrections of higher order in Za are indeed 
necessary to account for. 
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VI. FIGURE CAPTIONS 

1. Diagram for pair production with capture in relativistic heavy ion colliders. 

2. Energy spectrum of the emitted positron in pair production with capture for 100 GeV/nucleon Pb + Pb collisions 
at RHIC in the frame of reference of the nucleus where the electron is captured. The differential cross section is 
multiplied by the electron rest mass. The cross section is given in barns. The solid line is the exact calculation, 
while the dashed line is the analytical approximation, eq. 

3. Total cross sections for pair production at RHIC as a function of the nuclear charge of the ions. The solid line 
is the exact calculation, and the dashed line is obtained from eq. ^ given in ref. 
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4. Cross sections for pair production with electron capture for Pb^'^'^ at 33 TeV impinging on several targets with 
charge Z. The recent experimental data of ref. is compared to eq. ^ (dashed line) and to eq. |3| (solid 
line). 
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